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Abatract 

A computer program to calcula'e tranalimt and 
ateady state U-mporaturcs, pressures, and coolant Hows 
In a cooled turbine blade or vane with an impingement 
insert is described. Input to the progrr.ni includes a 
descriplliai of the blade geometry, coolant supply condi- 
tiona, outside thermal boundary conditions and wheel 
speed. Coolant* side heat transfer coefficients are cal- 
culated Internally in the program, sdth the user specify- 
ing the mode of heat transfer at each Internal flow ata- 
tlon. Program output includes the temperature at each 
node, the coolant pressures and flow rates, and the in- 
side heat transfer coefficients. A sample problem is 
discussed. 

Nomenclature 

b isiuivalent impingement slot width 

CI-C7 user supplied constants in general linpinge- 
mimt correlatkNi 

Dj impingement hole diameter 

pin fin diameter 
P 

DI-D8 user supplied constants in general leading 
edge Impingement correlation 

h^ coolant side heat transfer coefficient 

h|^ hot gas side heat transfer coefficient 

L^ pin fin length 

I surface half length over which leading edge 

Impingement heat transfer is averaged 

M film cooling mass velocity ratio, c-xilant to 

free stream 

m exponent in cf|uatlon (1) 

NUq Nusselt number based on channel hydraulic 
cUumeter 

NU| Nusselt number based on Impingement jet 

diameter 

Pr ITandtl number 

Re|j Reynolds number based on channel hydraulic 

diameter 

Rcj Reynolds number based on hole diameter 

RC|| Reynolds number based on equivalent slot 

width 
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pin fin spacing 
Stanton number 

equivalent slot height for film cooling 
coolant crossflow velocity 
Impingement Jet velocity 
impingement hole spanwise spacing 
distance downsiresm of film cooling hole or 


impingement channel width 
I)' film cooling effectiveness 

outside gas vi. cosity 
Pj film cooling Jet viscosity 

crossflow gas density 
Pj Impingemimt Jet density 

V>p s >2 constants in equation (IJ 

Introduction 

As core turbine engine operating conditions become 
more severe. It bi-comcs more difficult to efftictively 
cool blades and vanes. Advanced transient thermal cil- 
cuLtt'.onal ti-chniqiK‘S are needed in order to design reli- 
able turbine bladt^s. There appears to be no generally 
available computer program that utilizes these advonci-d 
techniqucd in combining the required heat transfer and 
coolant flow distribution calculations. Hius, it was de- 
cided to create a computer program that would perform 
transient and steady state heat transfer and coolant flow 
analyses for a cooled blade, given the outside hot gas 
boundary conditions, the coolant Inlet pressure or flow 
rate, the geometry of the blade shell, and the cooling 
configuration. 

1310 resulting program, TACTI, can handle a tur- 
bine blade or vane which is i>quipped with a central cool- 
ant plenum insert, from which coolant air flows through 
holes to impinge on the inner surface of the blade shell, 
it Is assumed that the spent Impingement air then flows 
chordwise and Is dumped through a split or drilled trail- 
ing edge and/or through film cooling holes. Hic blade 
is modeled by dividing it by chordwise planes into slices, 
with e .ch slice broken up into chordwise calculational 
stations. Temperatures at each station are calculated 
for foul points through the wall and one in the coolant 
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chanMi. Included In this modal U tha capalditty to ana- 
lyse a blade with a ceramic thermal harriar coatln(. 

The TACTI program Is used at the NASA Lewis Re- 
search Center on an IBM T8ii/360-6T computer. The 
source program conalats of approximately MOO lines at 
code and the program requires about 00 000 words of 
atorage. Typical running limes for the program are in 
the neighborhood of 1,4 seconds of CPL' lime per calcu- 
Istlonal station (or a steady state nn, and nbout 0,4 sec- 
onds of CPr time per station per time step (or a tran- 
sient run. 

The TACTI program Is described In detail In Refs, 
t and 2, This paper presents an overview o( the program 
and a sample problem to lUuatrate the use and capabili- 
ties of Ihe program. 

Method of Analysis 
Blade Geometric Model 

The key to creating a usable computer program is 
to have as simple a geometric model as possible (or the 
system being analysed. In this program, the emphasis 
is on a blade or vane with a central coolant plenum and 
chordwiae flow of the coolant after impingement, there- 
fore, it was derided that the primary calculational direc- 
tion would also be chordwlse. The blade la divided Into 
layers which are bounded by chordwlse cuts through the 
blade as shown In KIg, 1, Each slice Is treated sepa- 
rately In the program, with radial heat conduction In the 
wall the only communication between layers. 

Figure 2 shows the details of the geometric model 
for a single slice o( a blade, showing the breakdown of 
the blade or vane Into calculational stations and nodes. 
Each calculational station Is broken down into five nodes, 
located at: (1| the outer surface, (2) the interface be- 
tween coating and blade metal, (3) a point midway 
through the wall metal, (4) the wall Inner surface, and 
(5) at a mid-coolant channel location. 

For input to the program, the following basic ele- 
ments of the geometry are needed (or each station: (1) 
the thickness at the wall coating and the wall metal and 
coolant channel width, (2) the chordwlse distance o( each 
node (com the adjacent lower numbered node, and (3) the 
radial span (or this slice. In addition, depending on the 
mode o( heat transfer specified, the user must supply 
diameter and spacing for impingement holes or the di- 
ameter and spacing for pin fins. The input Is described 
In detail In Ref. 1. 

Numerical Model 

The numerical solution for the temperatures 
throughout the blade involves writing a transient energy 


balance equtUon for each nude and forming a set of equa- 
tions to be solved for tha temperature distrlbutioa. 
Similarly, the coolant pressure dtstiibutlon Is deter- 
mined by writing the trinslent momentum equation for 
flow between adjacent fluid nodt s and solving the result- 
ing set of nc|uations for static pressures. 

The nodal energy balances are linearised one- 
dlmenslonal heal conduction equations at the outside nodi, 
the junction of cL^dding and metal node, and the Inside 
wall aurface node. At the mIdmeUi node, a linearised 
three-dimensional heat conducticsi equation Is used. In 
the coolant channel, energy and momentum equatiuns (or 
one- dimensional compressible flow Including friction and 
heat transfer are written (or Ihe elemental channel length 
between two coolant nodes. Hie equations used are pre- 
sented In Ref. 1. 

Heat Transfer Correlations 

Three different mooes of coolant side heat transfer 
are built Into the program. The user must indicate the 
mode to be used ,it each station. Built-in correlations 
are available for: (I) Impingement, including separate 
correlations for the stagnation point in the con?ave lead- 
ing edge and (or stations where crossflow Is present, (2) 
forced canviK:llon channel flow, and (3) forced convection 
over on equilateral triangular array of pin fins. In addi- 
tion, the program has two general correlations that may 
be used In place uf the specific impingement correlations 
by including Ute appropriate constants In the Input, 

For impingement cooling, the correlation due to 
Kercher and Tabakoff, Including the effects of cross- 
flow, is used. 

( 2 \0.001 

D^j 

where and m are functions of the ratio of jet 
spacing to jet diameter and y >2 Is a function of the ratio 
of coolant crossflow to jet flow. The valuer of the coef- 
ficients can be found in Ref. I. The alternative general 
correlation is of the form: 



( 2 ) 


where the constants Cl through C7 arc specified by the 
user In the input. A correlation specifically for lmplng<>- 
ment Into a concave surface Is used to calculate the cool- 
ant side heat transfer coefficient at the leading edge im- 
pingement stagnation point. The built-in correlation Is 


due to Metxaer, Ysnuiahlta, and JenklDa/ aitd ha* the 
form 

-0.27 /I V®‘ ** 


at • 0. aSi He. 


(3) 


where the Reynolda number la baaed on the equivalent 
alot width, b, and I la the aurfaro half lenuth over 
which the heat tranafer la averaaed. The alternate gen- 
eral correlation la of the form 



where the conatanta Dl through IV are apecifled b)' the 
uaer In the Input. 

The correlation uaed for turbulent forcod convection 
channel flow la: 


NUp - 0. 023 He^- » 333 ^ jj 

where the auuacrlpt D Indlcatea that the channel hydrau- 
lic diameter la to be uaed aa the reference length. 

The correlation for forced coovectioo In an equi- 
lateral triangular array of pin ftna la taken from 
Faulkner. ‘ and la: 

Nup - 0. 023 
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The program contalna the capability to handle local 
film cooling from a row of holea or a alot. However, 
due to a program requirement for continuoua rearward 
flow In the coolant channel, the uae of the film coolli.g 
option mual be limited. If local film cooling la Included, 
the uaer haa the option of specifying the oiitnldc heat 
trantfer coefficient for the film cooled cane directly, or 
specifying an unblown heat transfer coefficient and let- 
ting the program calculate an effectiveness. An effec- 
tiveness Is calculated from the correlation due to 

A 

Stollery and El- Ehwany: 



where a is aii equivalent slot height, x Is the distance 

from the row of holes, M is the mass velocity ratio, 

coolant to free strean RCj la the couluit Jot Reynolds 

nu.mber, and u, and u_ are rh> coolant and free 
J g 


stream vtscosltlea. 

The air properties needed la the various correla- 
tions are evaluated In the program at a local reference 
temperature. A table of air properties versus tempera- 
ture at 30 atmospheres pressure, from Ikiferl and 
Svehlo,^ Is Included In TACTl. 

l*rogram Deaciiptlon 

Figure 3 shows on overall flow chart for the TACTl 
program. Hie method of aolutlon can be shown by de- 
scribing what Is involved In doing a steady state calcula- 
tion for a multi- slice blade. For transient runs, a 
steady state solutioa must be done first to provide the 
initial stste for the transient. There are three basic 
nested calculatlonal loops that must converge for a solu- 
tion to be reacbeiL These are labeled as A, B, and C 
In Fig. 3, The Innermost loop. A, results In a stable 
set of temperatures at all nodes and pressures at the 
coolant channel nodes. The intermld'ate loop, B, re- 
sults In a stable coolant flow split between suction and 
pressure sides of the blade, and the outermost loop, C, 

Is an overall coolant mass balance birtween tlie coolant 
sigiply and the coolant discharge to the main stream. 

The program starts with the coolant aigg>ly pressure 
and total coolant flow fixed. The Impingement flow Is 
Initially assumed to split uniformly at the leading edge 
stagnation station, station 1. All coolant flows for the 
slice laider consideration are calculated first, baaed on 
the latest pressure distribution. The temperatures at 
each node are then calculated by solving slmtiltam’ously 
the energy equations presented in Ref. I, and the pres- 
sures at all coolant nodes arc calculated by solving si- 
multaneously the momentum cKiuatloni presented in 
Ref. 1, This cycle of calculating coolant flows, all 
temperatures, and coolant channel pressures Is repeat- 
ed until the pressures converge. The flow split between 
suction and pressure side coolant channels Is then 
checkfKl by comparing the pressures at the entrance to 
the trailing edge region. If ttu>y do not match, the Im- 
pingement flow split at the leading edge Is adjusted and 
the Inner loop calculations are repeated. Once the 
proper flow split Is achieved, the program moves up the 
blade to the next slice and repeats the above sequence. 
After all N slices have converged, Uu total coolant 
flow used Is compared to the inlet value. If ther<' Is on 
imbalance, either the inlet flow or supply pressure Is 
adjusted, tk-pcndlng on 'vhlch was specified In Ihe input, 
s-.d the calculations start over again, tinre the u\eral! 

I -silont balance Is satisfied, the steaily state solution Is 
complete aiid the transient calculations begin. During 
a trar.sicmt, loop B Is bypassed, as the coolant flow 
spilt Is primarily a function of blade geometry. 
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A concepCual desipi ot an Implnfcpinntt cooled blade 
(or an advanced hiKh preaaure turbine la uaed to demon- 
strate the program. Two casea were nai, a rapid cnulno 
accel<<ratlan transient, and a almilar onipne deceleration. 
The transient calculations were carried (Hit to S seconds, 
usinx 0. 25 second time steps. Hiauilnx times on the 
TB8/MO-67 computer were show 3200 seconds at CPI' 
time (or each 5 second transient. 

Ihc blade considered has a span at 3. 81 cm and is 
divided Into Uirec slices. FiRure 4 is a mid-span cross 
soc..on view at (he blade showinR the locations at the 
computational stations and impingement holes. Figure 5 
shows the time variations assumed (or the transiimt In- 
put variables and Fig, 6 Is a plot o( the hot gas side heat 
tninatcr coe((iclent distribution (or these esses. The 
heat trans(er coe((iclents were calculated (rom a cylin- 
der leading edge correlation and the STANS boundary 
layer computer program at Re(. 8. 

Figures 7 to 9 are representative output (rom the 
TACTl computations. Figure 7 shows h^, the coolant 
side heat trans(cr coefficient distribution calculou-d by 
the program (or the high speed steady state starting 
condition. Note the very high host trana(er coetdeients 
around the (orward stagnation point due to pure Impinge- 
r.ent into a concave su'^ace. Further down the blade the 
h|, values drop atf as the effect at coolant crosadow Is 
(cit. The sharp dip In the curve (or both suction and 
pressure sides occurs at the end at the impingement 
insert where there Is a reduction In coolant velocity. 

Figure 8 shows the high speed steady state tempera- 
ture distribution around the blade (or a slice at the blade 
hub region, a slice at the blade mid- span region, and a 
slice at the blade tip region. The temperatures plotted 
'n Fig. 8 are mid-wall temperatures. Although the hub 
and tip regions have the same hot gas conditions and 
cooling condgurations (or this pisiblem, the tip region 
runs cooler due to the coolant having a higher pressure 
In this region. 

Figure 9 presents the transient behavior ot some o( 
the temperatures. The three curves are (or the outside 
sur(acc, mid- span hot spot, the outside surface, mid- 
span cold spot, and the overall bulk metal temperature. 
ITie symbcls on the curves Indicate the location on the 
blade of the hot and cold spr*s. Note that the locations 
of the mid- span hot spot and cold spot change during the 
transient. 
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Figure 2 . ~ Blade geometric model. 
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Fiqurfe 5. - Transient input for sample problem. 
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Figure 8. - Mid-wall steady state temperature distribution. 




Figure 9. * Transient behavior of blade temperatures. 
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